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Abstract 

We investigate the Painleve analysis for a (2 + 1) dimensional Camassa-Holm 
equation. Our results show that it admits only weak Painleve expansions. This then 
confirms the limitations of the Painleve test as a test for complete integrability when 
applied to non-semilinear partial differential equations. 
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1 Introduction 



Recently there has been a great deal of interest in the study of the integrability 
properties of the Camassa-Holm (CH) equation, 

u t + 2ku x - u xxt + 3uu x = 2u x u xx + uu xxx , (1) 

where subscripts denote partial derivatives and k is a constant [1]. Equation 
(1) was first derived by Fokas and Fuchssteiner using the method of recursion 
operators [2]. However, it became more widely known when it was derived from 
physical considerations as a water wave equation by Camassa and Holm, by 
using asymptotic expansions directly in the Hamiltonian for Euler's equations 
governing inviscid incompressible flow in the shallow water regime [1]. They 
showed that equation (1) is completely integrable for arbitrary values of k, 
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and that for k = it admits a special kind of travelling wave solution, the 
peakon, of the form ce~^ x ~ ct ^. Further, they showed that it possesses a Lax pair 
and bi-Hamiltonian structure. Subsequently, several works have been devoted 
to studying the underlying mathematical and physical properties of equation 
(1) [3,4,5,6,7,8,9,10,11,12,13]. 

Interestingly, several higher dimensional generalizations of the CH equation 
have also been proposed. For example, Kraenkel and Zenchuk constructed the 
following (2 + 1) dimensional integrable generalization of the CH equation [14], 



m t — um x — 2mu x = 0, 
mu y + u x - u xx - 2p x = 0, 

m 2 py + mp x + mp xx - pm x - m x p x = 0, (2) 

where m, u and p are functions of (x,y,t). In [14], equation (2) was derived 
from a system of (2 + 1) dimensional shallow water equations by using a 
multiscale decomposition, and was also shown to have a Lax pair. The spectral 
problem for this equation was solved by Zenchuk [15]. Kraenkel et al. [16] 
studied the invariance properties of (2) through Lie symmetry analysis and 
explored its similarity reductions, giving Lax pairs for all the resulting (1 + 1) 
dimensional partial differential equations (PDEs) and thus establishing their 
integrability. Recently, Johnson has derived another (2 + 1) dimensional CH 
equation, which has a structure similar to that of the Kadomtsev-Petviashvili 
equation, again within the context of equations governing water waves [17]. 
The integrability of this equation is discussed in [18] . 

Even though many interesting mathematical properties have been studied for 
the (2 + 1) dimensional CH equation (2), an important question remains unan- 
swered, that is, whether or not it passes the Painleve (P-) test [19,20]. Indeed, 
certain difficulties arise when applying the P-test to equation (2). These diffi- 
culties are similar to those that arise for the class of equations studied in [4], 
and which also arise for the special case of the CH equation itself. However, 
in [4] it was shown how these difficulties can be overcome by including an 
extra lower order term in the Painleve expansion; in particular, it was shown 
that the CH equation admits only weak Painleve expansions. We find that a 
similar modification of the leading order analysis allows us to apply the P- 
test to equation (2). We find that this equation admits only weak Painleve 
expansions. 
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2 Painleve analysis 



In this section we apply the P-test [19,20] to the (2 + 1) dimensional CH 
equation (2). The P-test consists essentially of three steps: (i) determination 
of leading-order behaviours, (ii) identifying the resonances, and (hi) check- 
ing the corresponding compatibility conditions (in order to avoid logarithmic 
branching) . 

2.1 (i) Leading order analysis 

As a first step we consider seeking local Laurent expansions about a nonchar- 
acteristic movable singular manifold <p = 0. We find that the leading order 
exponent for u cannot be negative, and so assume, following [4], that the 
leading orders of the solutions of equation (2) are 

u ~ cipt + Uo<fi a , manioc/) 13 , p~p </> 7 , (3) 

where c is a constant to be determined. We note that here we are using 
Kruskal's ansatz [21], that is, we are taking (j) = x + ip(y,t) and all coeffi- 
cients in our Laurent expansions (for example, u , m and p above) to be 
functions of (y, t) only. 

Balancing the most dominant terms, we find that c = 1, and are led to the 
following two cases: 

Case 1 a — |, f3 — —1, 7 = —\ with m = ^r, po = —^f, u o arbitrary. 
Case 2 a — |, f3 — —1, 7 = \ with mo = — ^r, «o and po arbitrary. 

2.2 (ii) Resonances 

The next step in the P-test is to find the resonances, that is, the powers at 
which arbitrary functions appear in the series. As we have encountered two 
branches in the leading order analysis, we present the calculations for each 
case separately. 

Case 1 

Substituting the expressions 
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m = mo^" 1 + m r (f) r ~ 1 , 

p = + Pr <f> r -*, (4) 

into the dominant terms tells us that the recursion relation for the coefficients 
of our Laurent expansions is of the form 



/ 


r 


n \ 

ru (J 








_i(2r + l)(r-l) 


lu if)y 1 - 2r 




= 9r 


V 





l/o 1 \ r(2r-l) 

s< 2r V, / 







where as usual f r , g r and h r are functions of ip, previous coefficients, and 
derivatives thereof. We thus obtain the resonances as r = — 1,0, J, The 
resonance r = — 1 is associated with the arbitrariness of the function if;; r = 
corresponds to the function u being arbitrary. The double resonance r = \ 
indicates that we must modify our Laurent series in order to include half- 
integer powers of 4>, thus obtaining a Puiseux series, and that two of the 
coefficients (mi , mi , pi ) therein will be arbitrary. Finally the resonance r = 1 
indicates that one of the coefficients (ui,mi,pi) will be arbitrary. 



Case 2 



In a similar way, substituting the expressions 



m = m o _1 + m r (f) r ~ l , 

P=P0<P^ +Pr4> r+ K 



(6) 



into the dominant terms of (2) tells us that the form of the recursion relation 
is 
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(7) 



where the functions f r , g r and h r are functions of ip, previous coefficients, and 
derivatives thereof. We thus obtain the resonances r = -1,-|,0,0,|. The 
double resonance r = corresponds to the arbitrariness of the functions u 
and p . The resonance r = | indicates once again that we need to modify our 
Laurent series in order to include half-integer powers of 4>, and that one of 



4 



the coefficients (ui,mi,pi) in the resulting series will be arbitrary. We might 
expect that the resonance r = — | means that we need to use the perturbative 
Painleve test. However, we will see that this resonance, on this occasion, can 
be accommodated much more easily. 

We note that, since we have non-integer resonances (in fact, in each of the 
above cases), equation (2) does not pass the Painleve test. It remains to check 
whether or not it admits weak Painleve — or Puiseux — expansions (see [20], 
also [4,22]). 



2.3 (Hi) Checking compatibility conditions 



We now check the compatibility conditions, using the following expansions. 



Case 1 



We substitute 



u = if) t + 2 ^2u 1 <f)2 i m = (j) 1 ^2m L (f)2, p = (p 2 ^2Pi<P 2 , (8) 
j=0 2 j=0 2 j=0 2 

with = x + 4)(y,t) and coefficients Ui = Ui(y,t), rrij = rrij(y,t) and p k = 
p k (y,t), and with m = p = — ^ and u arbitrary, into equation (2). 

At the point of determining the coefficients (mi, mi, pi), we find that pi is 

2 2 2 2 

arbitrary, and that one of ui and mi is also arbitrary, with both compatibility 
conditions satisfied. We choose to leave mi arbitrary, solving for ui as 

ui = — — Unmiih v . (9) 

At the point of determining the coefficients (ui,mi,pi), we find that any 
one of these three can be left arbitrary with satisfied compatibility condition. 
Thus for example we may choose m 1 to be arbitrary, and solve for Ui and p\. 
All subsequent coefficients are then determined in terms of the five arbitrary 
functions (if), Uo, mi ,pi , mi). Since all compatibility conditions are satisfied, 
we see that in this case, equation (2) admits a weak Painleve expansion. 



Case 2 



In this case, instead of the series that would usually be chosen to check com- 
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patibility conditions, we substitute the series 

oo oo ^ oo 

u = tp t + 02 u i4> 2 , m = <p^ 1 ^2m l (p2 ; p = p _i + (f)2 Ypi^, (10) 
j=0 2 j=0 2 j=0 2 

with m = — 2^-, and u and p arbitrary, into equation (2). The difference 
with the usual series that would be taken lies in the inclusion of the extra 
term p_i. It is the inclusion of this extra term that allows us to accommodate 

2 

the resonance r = — |, since we find that p_i is left arbitrary. 

At the point of determining the coefficients (ui,mi,pi), we find that any 
one of these three can be left arbitrary with satisfied compatibility condition. 
Thus for example we may choose mi to be arbitrary, and solve for ui and pi. 
All subsequent coefficients are then determined in terms of the five arbitrary 
functions (ip,p_i,uo,po,mi). Since all compatibility conditions are satisfied, 
we see that also in this case, equation (2) admits a weak Painleve expansion. 



3 Conclusions 

In this letter we have carried out a detailed investigation of the Painleve 
analysis of the completely integrable (2 + 1) dimensional generalization of the 
CH equation, equation (2). We have used a modified version of the leading 
order analysis, as adopted in [4] . In this way, we have shown that the equation 
under study admits only weak Painleve expansions. This example therefore 
confirms the limitation of the Painleve test as a test for complete integrability 
when applied to non-semilinear PDEs. 



Acknowledgements 

The work of MS forms part of a Department of Science and Technology, Gov- 
ernment of India sponsored research project. The work of PRG and AP is 
supported in part by the DGESYC under contract BFM2002-02609, and that 
of AP by the Junta de Castilla y Leon under contract SAO 11/04. PRG cur- 
rently holds a Ramon y Cajal research fellowship awarded by the Ministry of 
Science and Technology of Spain, which support is gratefully acknowledged. 
This work was carried out during the visit of MS to the University of Sala- 
manca, funded under contract BFM2002-02609, in June 2003. 



6 



References 

[I] R. Camassa and D. D. Holm, Phys. Rev. Letts. 71, 1661, (1993); 

R. Camassa, D. D. Holm and J. M. Hyman, Adv. Appl. Mech. 31, 1, (1994). 

[2] B. Fuchssteiner and A. S. Fokas, Physica D4, 47, (1981). 

[3] F. Cooper and H. Shepard, Phys. Lett. A 194, 246, (1994). 

[4] C. Gilson and A. Pickering, J. Phys. A: Math. Gen. 28, 2871, (1995). 

[5] B. Fuchssteiner, Physica D 95, 229, (1996). 

[6] M. Senthil Velan and M. Lakshmanan, Int. J. Nonlinear Mech. 31, 339, (1996). 

[7] J. Schiff, Physica D 121, 24, (1998). 

[8] M. Fisher and J. Schiff, Phys. Lett. A 259, 371, (1999). 

[9] A. Constantin, R. Soc. Lond. Proc. Ser. A Math. Phys. Eng. Sci 457, 953, 
(2001). 

[10] Z. Liu and J. Qian, Int. J. Bifurcation and Chaos 11, 781, (2001). 

[II] A. Pickering, Theor. Math. Phys. 135, 638, (2003). 

[12] R. Camassa, Discrete and Continuous Dynamical Systems 3, 115, (2003). 

[13] H. R. Dullin, G. A. Gottwald and D. D. Holm, Fluid Dynamics Research 33, 
73, (2003). 

[14] R. A. Kraenkel and A. I. Zenchuk, Phys. Lett. A 260, 218, (1999). 
[15] A. I. Zenchuk, Physica D 152-153, 178, (2001). 

[16] R. A. Kraenkel, M. Senthilvelan and A. I. Zenchuk, Phys. Lett. A 273, 183, 
(2000). 

[17] R. S. Johnson, J. Fluid Mech. 455, 63, (2002). 

[18] P. R. Gordoa, A. Pickering and M. Senthilvelan, Z. Naturforsch. A 59, 640, 
(2004). 

[19] J. Weiss, M. Tabor and G. Carnevale, J. Math. Phys. 24, 522, (1983). 
[20] J. Weiss, J. Math. Phys. 24, 1405, (1983). 

[21] M. Jimbo, M. D. Kruskal and T. Miwa, Phys. Lett. A 92, 59, (1982). 
[22] A. Pickering, Phys. Lett. A 221, 174, (1996). 



7 



